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Abstract. On the basis of integral representations of Poisson and binomial distri- 
bution functions via complete and incomplete Euler V- and B-functions, we introduce 
and discuss continuous counterparts of the Poisson and binomial distributions. The 
former turns out to be closely related to classical Volterra functions as well. Under 
usual condition Np — > A, we also prove that the sequence of continuous binomial 
distributions converges weakly to the continuous Poisson one. At the end, we discuss 
a relationship between the continuous Poisson distribution and the r-process. 



1. Introduction 

In various applied research papers, many authors extensively use what they call a "con- 
tinuous Poisson distribution" and a "continuous binomial distribution" , providing these 
terms with very different, not always correct meanings. For example, by the "continu- 
ous Poisson distribution" , some authors ( [8] , [12] ) understand an absolutely continuous 
distribution with the density of the form 

e~ x \ x 

f x (x) = C\— — — -, x > 0, 
L(x + 1) 

where C\ is a normalizing constant. Some others use to this end the T-distribution (pQ, 
[6]) or even simply the exponential distribution (|14j). 

From the strictly formal point of view, the above distributions can not be regarded 
as genuine continuous analogues of the classical Poisson law, since these distributions 
have probabilistically little in common with this law. In this paper, we make an attempt 
to introduce continuous counterparts of Poisson and binomial distributions in a more 
natural way, on the basis of integral representations of Poisson and binomial distribution 
functions via complete and incomplete Euler T- and B-functions, and examine some 
properties of these counterparts. 

Note that these distributions appeared in various forms (although without any detailed 
exposition) in the papers of G. Marsaglia [10] , H. Stern [12], as well as in the solution 
of problem 3.4.1.22 in the book of D. Knuth [9]. In paper [10] and book [9] they were 
employed for the computer simulation of the classical Poisson and binomial distributions. 
In report [12], it is noted that the continuous Poisson distribution is closely related to 
the distribution of the time when r-process jumps over a fixed level. This relation will 
be discussed in detail at the end of the paper. 

In the present work, we proceed with the research started in paper [7]. 
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2. Notation and preliminaries 
Recall the standard notation which will be used in what follows: for x G M set 
\x\ = f max{fc G Z: k < x}, \x\ = f minjfc G Z: k > x}. 
Besides, we will need incomplete Euler T- and B-functions: 

/•OO 

T(x, X) = e - *^" 1 dt, x > 0, A > 0, (incomplete T-function) 

B{x, y,p) d = / t x - x (l - tf- 1 dt, x,y>0,0<p<l. (incomplete B-function) 
J P 

In particular, we have for the usual (complete) T- and B-functions: 
r(i) = T(x, 0), B(x, y) = B(x, y, 0), x, y > 0. 

Denote by 7Ta, A > 0, the classical Poisson distribution, i.e. a discrete probability 
measure on the set N U {0} given by n\({k}) = e~ A ^-. Also denote by Pn, p , N G N, 
p G [0, 1], the classical binomial distribution, i.e. a discrete probability measure on the 
set {0, . . . , N} given by N , p ({k}) = ( N k )p k (l - p) N ~ k . 

It is well-known that the distribution function F\ of the Poisson measure -k\ admits 
a representation in terms of functions T(-, A) and T(-): 



F x (x) d =7rx((-oo,x)) = \ V, _ M _ lA » = <j r(w,x) (1) 

A similar representation via B-functions holds true for the distribution function -Fjv.p of 
the binomial measure /3n, p - 

!0, z<0, 
1, x > N, 

'0, x < 0, ( 2 ) 

BtM.AT+l-M.p) 0<x<Ar; 
B([x],iV+l-M) 

1, x > N. 

Representations ((T|) and ([2]) can be attributed to a sort of "probabilistic folklore" 
(see e.g. where they are given as problems). The proofs are a simple exercise in 

integration by parts. 



'0, a;<0, 



3. Continuous analogues of Poisson and binomial distributions 

Taking representations (JXJ) and © into account, the following definitions appear to 
be quite natural. 

Definition 3.1. By continuous Poisson distribution with parameter A > we will mean 
the probability measure fr\ supported by [0, oo) with distribution function of the form 



Fx(x) = rLjo _ " ' (3) 
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Definition 3.2. By continuous binomial distribution with parameters N > 0, p G (0, 1) 
we will mean the probability measure /3at iP supported by [0, N + 1] with distribution 
function of the form 



Fn, p (x) = < 



0, x < 0, 

B(x,N+l-x,p) Q < N 1 ,^ 

1, a:>iV + l. 

Note that in Dehnition 13.21 parameter TV may take non-integer values; we continue to 
denote it with N only by tradition. 

We should hrst check whether the distributions introduced by Definitions 13.11 and 13.21 
are well defined. To this end, we show that the following assertions hold true: 



(i 
(ii 
(hi 
(iv 

(v, 
(vi 



lim^+o Ej^y = 0; 
iiiAi x— f+oo r(x) ' 

the function ^jfej^ increases in x on the interval [0, oo); 

B(x,N+l-x,p) _ n . 
,. B(x,N+l-x,p) _ -, 

iim^^jv+i-o B (x,n+i-x) - l > 

the function B q X ^ N ^i^^ increases in x on the segment [0, N + 1]. 



Assertions (i) and (iv) hold, since the numerators are bounded as x —> +0 while the 
denominators increase to infinity. Similarly, B(.t, iV + 1 — x) — B(x, N+l — x,p) is bounded 
as x — > N + 1 — while B(ir, N + 1 — x) increases to infinity, which yields (v). (ii) follows 
from the fact that T(x) — T(x, A) has an asymptotic order x~ x X x as x — >• +oo, in other 
words is o\T(x)). 

Finally, assertions (hi) and (vi) result from the following relations: 



oo />oo g 



h(x)^F^(x)=r- 2 (x)[ I I e-^Hsty-Hn^dsdt 



A 

T- 2 {x)(^J J e-^+Q (st)*- 1 In ^ ds dt) >0, x > 0, 



(5) 



B- 2 (.x,7V + l-x)( / / (st^- 1 (I - s) N - x (l - t) N ~ x In y *? ds dt 
lp Jo *(1 - s) 

= B- 2 (a;,7V + l-a;)( / Astf'^l - s) N ' x {\ ~ t) N ~ x In ^ ~ datft) > 0, 
P 7o H 1 - s ) J 

< x < N + 1. 



(6) 



Note that we changed oo for A and 1 for p in the upper limits of the inner integrals due 
to antisymmetry of the integrands with respect to the argument pair (s,t). 

The result can be summarized in the next statement. Here R + stands for [0, oo). 

Theorem 3.1. The distribution functions F\ and Fn, p given by (0) and define 
absolutely continuous (with respect to the Lebesgue measure) probability distributions on 
the measurable space 23 (K + )) . Their densities are of form J2J) and fcfij). 
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4. VOLTERRA FUNCTIONS AND THE MOMENTS 
OF THE CONTINUOUS POISSON DISTRIBUTION 

In this section, we study the moments of the continuous Poisson distribution. They 
turn out to be closely related to classical Volterra functions, which play a role in integral 
transforms. 

We first give the corresponding definitions (see e.g. [5] or [3]). In the most general 
form, the Volterra /i-function is defined by 

f°° t x+a x^ 
fj,(t,a,P)= , , — — tt dx, Re/3>-l. 

Jo r (P + l)T(x + a + 1) 

In particular, by setting a — (3 — we come to a simpler and more known Volterra 
i/-function: v(t) = f Q r (* +1 ) dx. 

For a fixed k G N, consider the fc-order moment function of the continuous Poisson 
distribution: 

dcf / ~ 



r oa x k-x (t(x) — T(x, A)) 



TOfc(A) = / x k Tix(dx), A > 0. 

Then we have 

/>oo 

m k (X) = k -F\(z))dx = k I v ~ ^ . " y ~~'"" dx. 

Jo Jo r(z) 

By the definition of complete and incomplete T-functions and Fubini theorem 

r x ( f°° x k ~ x t x ~ x \ f x 

m h (X) = k e~* / -—dx)dt = k\ e'* fj,(t, -1, k - 1) dt. (7) 

Jo \Jo r (^) / Jo 

The moments of the continuous Poisson distribution look particularly nice in terms of 
Laplace transform in A. Again by Fubini theorem we have 



4 / \ def 

m k {s) = 



/ e- sX m k (X) dX = - e~* (a+1 W*, -1, k-1) dt. (8) 
Jo s Jo 

The Laplace transform fi of the Volterra function /i is of the form 

A(s,Q,/3) = ga +i (ms)/3 +i . Res>l. 

(See e.g. [3]; anyway, this can be deduced by a simple calculation.) Thus, flSJ) yields 

fc! 

S ln R (T+.s) 

The previous relation makes it possible to find the double Laplace transform (Laplace- 
Sticltjes with respect to measure and Laplace with respect to A) of the distribution family 
(7Ta,A > 0). Let £xj A > 0, denote a random variable distributed by the continuous 
Poisson law with parameter A. Then, for Re s > we have 

<p(p, s) d = f e- sX Ee- p ^ dX = [[ e -P^-^ f x (dx) dX = 

Jo J ./(0,ao)x(0,oo) 

i-P) k r „-.A„ V* (-P)*»M*) 1 ln(l + s) ( "" 



mfe(s) = fc -, Res>0. (9) 



EnF" / e - sA ^(A)dA = E 



fc=o ... „„ k=o ki sp + Hi + s y 

We now resume all the above in the next statement. 

Theorem 4.1. The moments of the continuous Poisson distribution tt\ are given by 
formula Ffy. In terms of Laplace transform they can be expressed as The double 

Laplace transform of the distribution family (tt\,X > 0) is of form 
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5. Convergence of continuous binomial distributions 

The classical Poisson theorem (sometimes also called the law of rare events) asserts 
that under the condition Np — > A the sequence of binomial distributions with parameters 
N,p weakly converges to the Poisson one with parameter A. In this section, we show 
that this remains true also for their continuous counterparts. 

Theorem 5.1. Let {fiN(k),p(k)i k E N) be a sequence of continuous binomial distribu- 
tions. If N(k) —> oo and p(k) — > as k — > oo in such a way that N(k)p(k) — > A then 

a w , ~ 

PN(k),p(k) — » 7TA- 

We break the proof up into three lemmas. 
Lemma 5.1. For arbitrary x G [0, AT] and y £ [0,oo) 

+ D) = rcx + im-x + i) '^ 1 -^"' (n) 

-A([y,2/ + l)) = T ^y- (12) 

Proof. Performing a simple integration by parts, we have 

B(x + 1,N - x,p) B(x,N+l-x,p) 



Pn, p ([x,x + 1)) = F N>p (x + 1) - F N , p (x) 



B(x + l,N-x) B(x,N + l-x) 
p x {l-p) N - x 



(N - x)B(x + 1, N - x) ' 

which yields (flTj) by the relationship between T- and B-functions. (fl"2f can be obtained 
in a similar manner. 



Lemma 5.2. Under the conditions of Theorem \5.1\ 

lim 0N(k), P (k){[x,x + !)) = ft* + 1)), x > 0. 



The proof follows from Lemma 15.11 and elementary computations with Stirling's for- 
mula and is thus omitted. 

Lemma 5.3. Let (Vk,k £ NU {oo}) be a sequence of probability measures on the mea- 
surable space (R + , Q3(R + )) . If for each x > 

lim V k ([x,x + 1)) =V 00 ([x,x + l)), (13) 
then the weak convergence of Vk to Voo holds. 

Lemma 15.31 actually implies that the intervals [x, x + 1), x > 0, form a convergence 
determining class (see e.g. [5], p. 18). 

Proof. Denote by Ff. the distribution function of the measure Vk- For a fixed e > 
let N e be such an integer that V 00 {[N e , oo)) < e. Then, for each k e NU {oo} and x > 

\N E -x] 

F k (x) = l-V k ([x,oo))=l- V k ([x + n-l,x + n))-V k ([x+\N e -x],(x>)), 

n=l 

and so 

\N e -x] \N e -x] 

Vao([x + n-l,x + n))- Fk([x + n- l,x + n)) + (14) 



< 



+^ BO ([a!+r^'e-*i»oo))+^fc([*+r jv =- aj i. 00 ))- 
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The absolute value on the right-hand side in (TT41 tends to zero as k — > oo due to 
condition (fT3|) . Since x + \N e —x] > N e , we have Poo ([a; + \N 5 — x] , oo)) < e. Moreover, 

limsupP fc ([a?+rJV e -a:l,oo)) < lim V k ([N s , oo)) = 1 - lim P fc ([0,iV e )) = 

= l-P oo ([0,JV e )) ='P 00 ([JV e ,c»)) <£• 

Letting now at first fc — > oo and then e — »• 0, we obtain by (fT4f that 

lim Ffc(a;) = ^(x), a; > 0, 

k— >oo 

which was to be proved. 

Proof of Theorem 15.11 follows immediately from Lemmas 15.21 and 15.31 



6. An application to the T-process 

In this section, we show that the continuous Poisson distribution appears as the dis- 
tribution of the (in a proper way normalized) time when T-process jumps over a fixed 
level. 

We start with an informal explanation of this phenomenon. Consider a sequence of 
independent identically exponentially distributed random variables {Ti,i £ N} and the 
sequence of their cumulative sums T = {T n ,n G N}: T n = X^ILi 7 *- Also set To = 0. 
It is well known that the renewal process [X(t),t > 0) constructed by the sequence T 
as X(t) — sup{n: T n < t} is the classical Poisson process. The question is whether this 
setting can be carried into continuous time in such a way that the continuous Poisson 
distribution occurs instead of the classical one. 

To answer this question we look upon the sequence T as a Levy process in discrete 
time, i.e. random sequence with independent homogeneous increments. We can carry this 
process into continuous time, since the distribution of unit increments (the exponential 
one) is infinitely divisible. What's more, it is a special case of the r-distribution. Thus 
we may regard the T-process as a natural continuous counterpart of the sequence T. 

Recall the corresponding definition (see e.g. [3J, P- 54). 

Definition 6.1. A Levy process \T(t),t > 0) is called the T-process with parameters 
a, /3 > if its transition density has a form 

Mx) = ^x at - l e-e*, x > 0. 

1 (at) 

Assertion 6.1. Let t c be a time when the process T jumps over a fixed level c. Then 
the random variable ar c is continuously Poisson distributed with parameter /3c. 

The proof follows from the fact that the distribution function F c of the random variable 
ar c coincides with that of the continuous Poisson distribution -kp c : 

F c (x) = P[r c < ^} = P{r(^) > C }^1° /*/*(«) du =^)J c ° C u x ~ l e-^ du = 
_ T(x,/3c) 



T(x) 



F Pc (x), x>0. 
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